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Abstract 



It is shown that the natural framework for the solutions of any Schrodinger 
equation whose spectrum has a continuous part is the Rigged Hilbert Space 
rather than just the Hilbert space. The difficulties of using only the Hilbert 
space to handle unbounded Schrodinger Hamiltonians whose spectrum has 
a continuous part are disclosed. Those difficulties are overcome by using 
an appropriate Rigged Hilbert Space (RHS). The RHS is able to associate 
an eigenket to each energy in the spectrum of the Hamiltonian, regardless 
of whether the energy belongs to the discrete or to the continuous part of 
the spectrum. The collection of eigenkets corresponding to both discrete 
and continuous spectra forms a basis system that can be used to expand 
any physical wave function. Thus the RHS treats discrete energies (discrete 
spectrum) and scattering energies (continuous spectrum) on the same footing. 
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Typeset using REVT^ 



1. INTRODUCTION 



Several authors have reahzed that the Hilbert space is not sufficient for the purposes 
of Quantum Mechanics, and that an extension of the Hilbert space to the Rigged Hilbert 
Space (RHS) is needed. The RHS was introduced in physics for the ffist time in the 1960s 
independently by J. P. Antoine [|T|, A. Bohm 0, and J. E. Roberts These authors realized 
that the RHS provides a rigorous mathematical rephrasing of Dirac's formalism. In essence, 
the Nuclear Spectral Theorem [^] (also known as Gelfand-Maurin Theorem) restates Dirac 
basis vector expansion along with the Dirac bras and kets within a mathematical theory. 
Later on, other authors such as A. Galindo and P. Pascual 0, and O. Melsheimer came 
to the same conclusion. 

Earlier attempts to go beyond the Hilbert space framework are reported in Refs. 
If two operators of the algebra of observables satisfy Heisenberg's commutation relation, 
at least one of them cannot be continuous (bounded) with respect to the Hilbert space 
topology. In Refs. , it is shown that there are subdomains of the Hilbert space that can 
be endowed with topologies that make those operators continuous. 

Since the early 1990s, the RHS has become a standard tool in many areas of theoretical 
physics, especially in those that deal with continuous and resonance spectra. For instance, 
the RHS has been used to treat the Lippmann-Schwinger equation P,[lO| , the Gamow vectors 



(cf. P,p!l|-p!^ and references therein), and certain generalized spectral decompositions of 
chaotic maps ||15|,|16|. 

The dynamical equation of Quantum Mechanics is the Schrodinger equation. Thus any 
attempt to show that the RHS contains the mathematical methods needed by Quantum Me- 
chanics should show that the natural framework for the solutions of the Schrodinger equation 
is the RHS. The objective of this paper is to show that the solutions of the Schrodinger equa- 
tion fall in a RHS. To illustrate this, we shall use the example of the square well-barrier 
potential Hamiltonian. We recall that none of Refs. |]TO|-|T^ took the Schrodinger equation 
as the dynamical equation. 

A RHS is a triplet of spaces 

$c7i:c$\ (1.1) 

where H is a Hilbert space, $ is a dense subspace of the Hilbert space, and is the dual 
space of i.e., is the set of antilinear functionals over the space $ . The space $ 
has a topology that is finer than the topology inherited from Ti. The domain 'D{H) of the 
Hamiltonian lies between $ and Ti, 

^CV{H)cn. (1.2) 

If a Schrodinger Hamiltonian is defined on the whole Hilbert space and has only discrete 
spectrum, then the Hilbert space is sufficient for the purposes of Quantum Mechanics. How- 
ever, if the Hamiltonian is not defined on the whole Hilbert space (e.g., it is unbounded) and 
its spectrum has a continuous part, then the mathematical methods of the Hilbert space are 
not sufficient, and an extension of those methods is needed. The RHS arises as the natural 
extension when dealing with unbounded operators that have continuous spectrum. 

The reason why the RHS is the extension of the Hilbert space that we need is twofold. 
On the one hand, one of the key assumptions of Quantum Mechanics is that the quantity 
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{if, Hip) 



(1.3) 



accounts for the expectation value of the measurement of the observable H in the state 
and that 

/\^H = ^{^,m^)-{^,H^f (1.4) 

accounts for the uncertainty of the measurement of the observable H in the state (p. The 
expectation value ( [1.3|) cannot be computed for every element of the Hilbert space 7i, but 
only for those if E Ti. that also belong to T>[H). Similarly, the uncertainty ( |1.4| ) cannot 
be computed for every element of H either. If we take as physical states those square 
integrable functions for which physical quantities such as the expectation value ( |1.3| ) and 



the uncertainty ( |1.4| ) can be computed, then not every square integrable function (i.e., every 
element of Ti) can represent a physical state. As we shall see, the natural space of physical 
wave functions is the space $ in (p..l|), because all physical quantities such as expectation 
values and uncertainties can be computed for its elements. 

On the other hand, in Quantum Mechanics it is assumed that for each energy in the 
spectrum of the Hamiltonian H, there corresponds a ket that is an eigenvector of H. If we 
denote the ket corresponding to an energy En in the discrete spectrum by \En) and denote 
the ket corresponding to an energy E in the continuous spectrum by \E), then we should 
have 

H\En) = En\E„) , (1.5a) 
H\E) = E\E). (1.5b) 

These eigenkets are normalized according to the following rule: 

(E„|E^)=5„^, {E\E') = 6{E - E') , (E|E„) = 0. (1.6) 

It is also assumed that those kets form a complete basis that can be used to expand any 
wave function ip, 



Y,\En){EnW)+ j dE\E){E\cp). (1.7) 



The solutions of Eq. (|1.5a|) are square normalizable, i.e., they lie in the Hilbert space. How- 
ever, the solutions of Eq. ( |1.5b| ) are not square normalizable, i.e., they lie outside the Hilbert 
space. Therefore, the Hilbert space is not large enough to contain the non-normalizable 
eigenkets that are associated to the energies in the continuous spectrum — a larger space 
than the Hilbert space is needed. This larger space, which contains the eigenkets \E), is the 
space $^ in ( |1 . 1|) . The action of the Hamiltonian, which is in principle defined only on the 
elements of $ (or on the elements of V{H)), can be extended to the elements \F) of by 
defining the following conjugate operator H^: 

i^lH^'lF) := {H^^\F) , V(^G$,|F)G$^ (1.8) 



The operator is a uniquely defined extension of H. Using definition ( [L.8| ), the RHS 
formalism restates Eq. (|1.5b|) as 
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{ip\H'' \E) = (H^iflE) = E{ip\E) , Vyp e * 



(1.9) 



When the arbitrary (p E ^ is omitted in this equation, we recover the formal equation ( |1.5bD . 



In this way, the Gelfand triplet ( |1.1D arises in a natural way. The Hilbert space Ti 
appears because the wave functions must be square normalizable. The subspace # is the 
set of physical wave functions, i.e., the set of square integrable functions for which any 
expectation value and any uncertainty can be computed. The dual space contains the 
eigenkets associated to energies in the continuous part of the spectrum of H. These eigenkets 
are defined as functionals over the space and they can be used to expand any G $ 
as in Eq. (1/7). The eigenequation ( 1.5b ) holds within the RHS in the sense of Eq. (1.9). 



Thus the Rigged Hilbert Space is the extension of the Hilbert space that is needed when 
the solutions of the Schrodinger equation are not normalizable. These non-normalizable 
solutions are treated as distributions in the dual space 

Unlike the Hilbert space, the RHS treats the discrete and continuous spectra on the same 
footing: for each energy of the spectrum of the Hamiltonian (either in the discrete or in the 
continuous part), there is an eigenket of the Hamiltonian (given by Eq. ( |1.5aD in the discrete 



case, and by Eq. (|1.5b|) in the continuous case) that belongs to a basis that expands any 



physical wave function as in Eq. ( |1.7| ). 

The main shortcoming of the RHS theory is that it does not provide a prescription to 
construct the space $ and the eigenvectors \E). The general statement of the Nuclear 
Spectral Theorem [^] just assures the existence of the generalized eigenvectors \E). In that 
theorem, the space $ is assumed to be given beforehand. In this paper, the RHS of the 
square well-barrier will be constructed explicitly along with the eigenkets \E) by applying the 
method proposed in Ref. This method has been successfully applied to a Hamiltonian 
with purely continuous spectrum: the square barrier Hamiltonian |]T^. Here we apply 
it to a Hamiltonian with both discrete and continuous spectrum: the square well-barrier 
Hamiltonian. 

The steps of our method can be summarized as follows: 

(i) Construction of the self-adjoint Hamiltonian from the formal Schrodinger differential 
operator. 

(ii) Construction of the Green function and therewith of the resolvent. 

(iii) Computation of the spectrum of the Hamiltonian. 

(iv) Construction of the direct integral decomposition. 

(v) Construction of the Rigged Hilbert Space. 

Steps (i)-(iv) will be addressed by using the Sturm-Liouville theory JTSf. Step (v) will be 



handled by the methods proposed in Ref. p. 

In Section we construct the self-adjoint Hamiltonian from the formal Schrodinger 
differential operator and from a domain of the Hilbert space. In Section |^, we construct 
the Green function and the resolvent of H. In Section ^, we compute the spectrum of the 
Hamiltonian. In Section |^, we construct the direct integral decomposition of the Hilbert 
space generated by the Hamiltonian. In Section |^, the RHS of the square well-barrier 
Hamiltonian is constructed. 
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2. FORMAL DIFFERENTIAL OPERATOR AND SELF-ADJOINT 

HAMILTONIAN 



The first step toward the construction of a RHS associated to a Hamiltonian is to define 
that Hamiltonian. In non-relativistic quantum mechanics, the Hamiltonian is defined by the 
formal differential operator of the time-independent Schrodinger equation and by a Hilbert 
space domain on which the Hamiltonian acts. 

The time-independent Schrodinger equation reads in the position representation as 

/\ + V{x)\ (x\E) = E(x\E) , (2.1) 

2m J 

where A is the three-dimensional Laplacian and 

-Vi <r < a 

V{x) = V{r) = { V2 a<r<b (2.2) 
6 < r < 00 

is the well-square potential, Vi and V2 being two positive real numbers. The potential V{x) 
is rotational invariant, and spherical coordinates x = (r, 9, 0) can be used to write Eq. ( p.l|) 

as 

r+ ^ }\' + V{r) (r, 6, <j>\E, /, h) = E{r, 6, <j>\E, /, h) . (2.3) 



2m r dr"^ 2mr'^ 
By splitting radial and angular dependences, 

(r, 9, /, h) ^ {r\E)i {6, h) ^ ^Xiir; E)Y^^i9, 0), (2.4) 

where Yij.^{6, 0) are the spherical harmonics, we obtain for the radial part 

+ ^ , ' + V{r) xi{r; E) = Exi{r; E) . (2.5) 



2m dr"^ 2mr^ / 
For simplicity we study only the case / = 0. Writing xi=oij'i E) = x(^; E), 

■^^2+Vir?jxir;E) = Exir;E). (2.6) 

Therefore, the formal Schrodinger differential operator for the zero angular momentum case 
is 

"^-J^l^r.^vC-)- (2.7) 

In order to construct the domain 'D{H) of the Hamiltonian, we need first the Hilbert 
space that contains it. Clearly, the Hilbert space that fits the differential operator ( ^.7[ ) is 
L^([0, 00), dr). Now, each element / in the domain 'D{H) must be square integrable (see 
condition (|2.9a|) below), the action of h on f must be well-defined (see condition (|2.9b|) 
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below), and the action of /i on / must remain square integrable (see condition (|2.9c| ) below). 
Furthermore, the Hamiltonian H induced by h and by T>{H) must be self-adjoint. It is 
well-known that all the possible self-adjoint operators associated to the formal differential 
operator ( p.7|) are determined by the following condition (cf. page 1306): 

/(O) + a/'(0) = 0, -oo<a<oo. (2.8) 

Among all the possibilities in Eq. the one that is used in physics is a = 0, i.e., /(O) = 0. 
All the above mentioned conditions, which determine the domain V{H), can be written as 

/(r) eL2([0,oo),rfr), (2.9a) 

f{r)eAC'[0,oo), (2.9b) 

ihf)ir)eL\[0,oo),dr), (2.9c) 

/(0) = 0, (2.9d) 

where AC^([0, oo)) denotes the space of all functions / which have a continuous derivative 
in [0, oo), and for which /' is not only continuous but also absolutely continuous over each 
compact subinterval of [0, oo). Thus the domain of our Hamiltonian is 

ViH) = {fe L\[0, oo), dr)\fe AC'[0, oo), hf G ^^([0, oo), rfr), /(O) = 0} . (2.10) 

It follows that the Hamiltonian defined by 

{Hf){r):={hf){r), f e V{H) (2.11) 

is a well-defined self-adjoint linear operator. 



3. RESOLVENT AND GREEN FUNCTION 

The second step toward the construction of the RHS is to compute the resolvent and 
the Green function of the (self-adjoint) Hamiltonian. The resolvent operator of H is defined 
as the inverse of the operator {E — H) for those complex energies E for which that inverse 
exists. If we write the resolvent operator as an integral operator, the kernel of that integral 
operator is precisely the Green function, 

/"OO 

[{E-H)-'f]ir)= dsGir,s;E)fis). (3.1) 
Jo 

The Green function G(r, s; E) can be computed using the prescription provided by the 
following theorem (cf. Theorem XIII. 3. 16 of Ref. []TB[ and also Ref. [IU|): 

Theorem 1 Let H be the self-adjoint operator ( |2.11|) derived from the real formal 
differential operator ( p.7[ ) by the imposition of the boundary condition ( |2.9d| ). Let Im(ii^) 7^ 
0. Then there is exactly one solution x(r; E) of [h — E)a = square-integrable at and 
satisfying the boundary condition ( |2.9d| ), and exactly one solution G(r; E) of {h — E)a = 
square-integrable at infinity. The resolvent {E — H)~^ is an integral operator whose kernel 
(j(r, s; E) is given by 
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2m xir;E)eis;E) ^ ^ ^ 



Gir, s;E) = {Z ' / : (3-2) 



2m x{s;E)e{r;E) 



where W{x, ©) is the Wronskian of x and 0, 

W{x, 0) = X©' - X'Q ■ (3.3) 



In order to compute G{r, s; E), we shall divide the complex energy plane in three different 
regions and obtain G{r, s; E) for each region separately. In our calculations, we shall use the 
following branch of the square root function: 

: e C I - vr < arg(E) < vr} i — ^ {E eC\ - 7r/2 < arg(E) < 7r/2} . (3.4) 



3.1. Region Rfi{E) < 0, Im.{E) 7^ 

We first apply Theorem 1 to the region of the complex energy plane where Re(£') < 
and lm{E) 7^ 0. According to that theorem, the Green function reads 

'2m/n? x{r;E)0{s;E) ^ ^ ^ 



Gir, s;E) = i v^^gP ^ ^ ^ ME) < , Im(E) ^ . (3.5) 

In this equation, the eigenfunction x(r; E) satisfies the time-independent Schrodinger equa- 
tion ( p.6| ) subject to the following boundary conditions: 

X(0;E) = 0, (3.6a) 

x(r;E)GAC2([0,oo)), (3.6b) 

x(r; E) is square integrable at , (3.6c) 

which yield 
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, U'^(E+Vi)r - U2^{E+Vi)r-. „ 

rgV 1' _ e V ^ < r < a 



X{r-,E) = \ J,(E)eV-1^(^-^-)- + J,(E)e-V"1^^^-^-)- a<r<6 (3-7) 
J3(^)e\/"^^" + J4(^)e^V^"^^" 6<r<oo. 

The expressions of the functions Ji{E)-Ji{E) can be found in Appendix 0. 

The eigenfunction Q{r\E) of Eq. ( ^.5] ) satisfies the time-independent Schrodinger equa- 
tion ( p.6| ) subject to the following boundary conditions: 

e(r;E) G AC2([0,oo)), (3.8a) 
9(r; £■) is square integrable at 00 , (3.8b) 

which yield 
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ME)e^-_^^' + ME)e-^~_^^- 0<r<a 
Qir-E)=\ l3(E)eV^-^(^-^^)'- + A(E)eV-1ff(^-^^)'- a<r<b (3-9) 
eVi^ b < r < oo . 

The expressions of the functions Ai{E)-Aa{E) can be found in Appendix 

3.2. Region Re(£^) > 0, \m.{E) > 
In this region of the complex plane, Theorem 1 leads to the following Green function: 

G(r, s-E)={ v^gr ^SS,.E) ^ ME) > , lm{E) > . (3.10) 



In this equation, the eigenfunction x(r; E) satisfies the differential equation ( |2.ti| ) subject to 
the boundary conditions (|3.6|), 



x{r]E) 



sin(^f^(E + \/i)r) 0<r<a 

J-i(E)eV¥(^'- + J,^E)e^^^^' a<r<b (S-H) 

^ J3(^)eV^'''' + J4(^)e"V^''' 6<r<oo. 

The functions Ji{E)-Ji{E) are determined by the boundary conditions (|3.6| ), and their 
expressions are listed in Appendix 0. 

The eigenfunction 0+(r; E) satisfies the differential equation ( ^iGf ) subject to the bound- 
ary conditions (p.8|). 



< r < a 



0+(r;i?) = { ^3+(E)_eV^(^-^^)^ + ^+(E)e-V^(^-^^)'- a < r < 6 (3-12) 
eVfi^ o<r<oo. 

The functions Af{E)-A^{E) are determined by the boundary conditions ( p.8|) , and their 
expressions are listed in Appendix 0. 

3.3. Region Re(^) > 0, Im(^) < 

If Re(-E') > 0, Im(i?) < and we use the branch of the square root function (p. 41), then 
the Green function has the form 

2m/fi2 X{r;E)e-{s;E) ^ ^ ^ 

G{r,s;E)=l xjs^f.U ^ Re{E) > , lm{E) < . (3.13) 
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The eigenfunction x{i^'-i E) satisfies the time-independent Schrodinger equation ( |2.6| ) subject 
to the boundary conditions (|3.6| ). Thus x{f]E) is given by Eq. ( 3.11|) — hence the same 



symbol as in the previous subsection — although now E lies in the fourth quadrant of the 
complex energy plane. 

The eigenfunction 9_(r; E) satisfies the differential equation ( p. 61) subject to the bound- 
ary conditions (|3.8|), 



A^iE)e'rf^^ + ^-(E)e-V^(^^ < r < a 
Q-{r;E)=^ ^3(E)eV^(^-^^)^ + A~(E)e-V^(^-^^)'- a<r<b (3-14) 
evA-^ < r < oo . 

The functions Ai {E)-A^ (E) are determined by the boundary conditions ( p.8| ), and their 
expressions are listed in Appendix 0. 



4. SPECTRUM OF THE HAMILTONIAN 

The third step toward the construction of the RHS is to obtain the (Hilbert space) 
spectrum Sp(if) of the self-adjoint Hamiltonian H. Since H is self-adjoint, the spectrum 
must be real. In order to elucidate which real numbers are in Sp(if) we make use of 
Theorem 3 below. Before stating Theorem 3, we need to state another theorem, which 
provides the unitary operator U that will be used to diagonalize H (cf. Theorem XIII. 5. 13 
of Ref. PI): 



Theorem 2 (Weyl-Kodaira) Let h be the formally self-adjoint differential operator 
defined on the interval [0, oo). Let H be the self-adjoint operator ( |2.11|) . Let A be an 
open interval of the real axis, and suppose that there is given a set {ai{r]E), a2{r;E)} 
of functions, defined and continuous on (0, oo) x A, such that for each fixed E in A, 
{(Ji(r; E), (T2{f", E)} forms a basis for the space of solutions of ha = Ea. Then there exists 
a positive 2x2 matrix measure {pij} defined on A, such that the limit 



{Uf)i{E) := lim lim 



Od- 



f{r)ai{r;E)dr 



(4.1) 



exists in the topology of L^(A, {pj^}) for each / in L^([0, oo), dr) and defines an isometric 
isomorphism U of E(A)L^( [0, oo), (ir) onto L^(A, {pjj}), E(A) being the spectral projection 
associated to A. 

The spectral measures {pij} are provided by the following theorem (cf. Theorem XIII. 5. 18 
of Ref. 



Theorem 3 (Titchmarsh-Kodaira) Let A be an open interval of the real axis, O 
be an open set in the complex plane containing A, and Re(if) := C — Sp{H). Let 
{ai{r;E), a2{r;E)} be a set of functions which form a basis for the solutions of the equa- 
tion ha = Ea, E & O, and which are continuous on (0, oo) x O and analytically dependent 
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on E for E in O. Suppose that the kernel G{r,s;E) for the resolvent {E — H) ^ has a 
representation 



\Y:lj=,0±(E)a,(T-,E)as(s-,E) r>s, 

for all E in Re(if) fl O, and that {pij} is a positive matrix measure on A associated with H 
as in Theorem 2. Then the functions Of^ are analytic in Re(if) fl O, and given any bounded 
open interval (£'1, E2) C A, we have for 1 < j < 2, 

p,,((Ei, E,)) = lim,^o lime^o+ /| " ^^e) - ^.^(^ + 

(4.3) 

= lim,^o lime^o+ lE!J[(^tjiE - le) - Of.iE + te)]dE . 



Using these spectral measures we can compute the spectrum of H. This spectrum is the 
subset of the real line on which the Green function fails to be analytic. This non-analyticity 
of G{r,s;E) will be built into the functions 9f^{E) that appear in Theorem 3. From the 
expression of the Green function computed in Section it is clear that the subsets {—00, 0) 
and (0, 00) should be studied separately. We will denote either of these subsets by A. 



4.1. Negative Energy Real Line: A = (— cxd,0) 



We first take A from Theorem 3 to be {—00^ 0). We choose a basis for the space of solu- 
tions of the equation ha = Ea that is continuous on (0, 00) x A and analytically dependent 
on E as 



Oi{r-,E) 




{E+V-i_)r 



{E-V2)r 



''-(E+Vx)r 



+ B2{E)e V 



-^Er 



<r < a 
a < r < b 
b < r < 00 . 



a2{r;E) = Q{r;E). 



(4.4a) 



(4.4b) 



The functions B1-B4 are such that o"i(r; E) and its derivative are continuous at r = a and 
at r = 6. Their expressions are listed in Eq. (|A8|) of Appendix 0. The function 0(r; E) is 
given by Eq. (|3.9[). Obviously, 

(4.5) 



X(r; E) = ME)a^{r- E) + ME)a2{r- E) 
which along with Eq. ( |3.5| ) leads to 



G{t,s-E) 



2m /h^ 



^-2m/h^ E 2 



a,{r-E) + ^^a2{r-E) 



a2{s;E) 



ME) 

r <s, Re{E) < , lm{E) ^ . 



(4.6) 



Since 
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a2{s;E) = a2{s\E) 



(4.7) 



we can write Eq. 
G{r,s-E) = - 

By comparing Eqs. 



as 



2m/ 



ai{r- E)a2{s- E) 



^ , ME) 



a2{r; E)a2{s; E) 



JsiE) 

r < s, Re{E) < , lm{E) ^ . 



(4.8) 



and ( [4.8|) we see that 

2m/h^ 1 



^-2min? E 2 
2m./R^ 1 J4(g) 



Re(E) < , lm{E) ^ . 



(4.9) 



The functions 0^j{E) are analytic in a neighborhood of A = (— cx3,0) except at the energies 
for which JziE) vanishes. In this case, the function 6*22 (-^) ^ pole. Since 



(4.10) 



where J+{k) is the Jost function the poles of 622{E) correspond to the zeros of the Jost 
function on the negative real axis of the energy plane, or on the positive imaginary fc-axis. 
Those bound states have been calculated by Hogreve pO|] for different values of Vi, V2 in 
Eq. (|2.2| ). In general, there is a finite number of bound states whose energies -Ei, • • • , E]\[ are 
in the interval (—1^1,0). The corresponding normalized wave functions read |]18 



0„(r) = A^„e(r;EO, n=l,...,N, 



(4.11) 



where the square of Nn is given by the residue of the spectral function ^22 (-^) the bound 
state energy 



res 



[^2^2 (^)] 



E—En 



res 



2m/h'^ 1 Ji{E) 



-2m/h? E'2 J^{E) 



(4.12) 



E—En 



Therefore Ei, . . . , E^ are the only negative energies that belong to the spectrum of our 
Hamiltonian. 

It is worthwhile noting that the normalization ([4.12| ) provided by the residue of the spec- 



tral measure 622{.E) is the same (as it should!) as the normalization obtained in Refs. ||2T|-|23 
To see this, we denote the momentum of those bound states by kn 
n = 1 N. Then we have that 



i kr, 



Q{r;k^) 



and that 



res 



2m/h^ 1 Ji{E) 



-2m/h^E'2j^(E) 



E—En 



e(r; -ikn) = 6+(r; kn) , 

2m/ IJJE) 
^j2m/h?E'2J+{E) 



- ikfii 
(4.13) 



E—En 



^^es[S{k)]k^kn ' 

(4.14) 
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where S{k) is the S'-matrix in the /c-momentum representation and 

J4E) = 2tJ3{E). 



(4.15) 



Since the eigenfunction (pnir) of Eq. ( [4 .111 ) is normahzed to one, from Eqs. ( [4.12| )-( p:.14| ) it 
follows that 



—I 



res [S{k)\ 



(4.16) 



which is the normalization that appears in Refs. [pi| -p3|. (The normalization rule of 
Eq. (4.16) has been generalized to Gamow vectors in Refs. pl|-p6[].) 



4.2. Positive Energy Real Line: A = (0, oo) 



Now we study the case A = (0, oo). In order to be able to apply Theorem 3, we choose 
the following basis for the space of solutions of ha = Ea that is continuous on (0, oo) x A 
and analytically dependent on E: 



ai{r;E) = xir;E) , 



a2(r;E) 



cos 



2m , 



(E + V^r) 



C3(^)eVf?^^ + C4(^)e" 



<r < a 

a < r < b 
b < r < oo . 



(4.17a) 
(4.17b) 



The functions C1-C4, whose expressions are given by Eq. (|A9|) of Appendix 0, are such that 
(T2 and its derivative are continuous at r = a and at r = 6. The eigenfunction x(r; E) is 
given by Eq. ( |3.11|) . 

Eqs. ( P^ , (^J§ and {^^J7^ lead to 



e+(r;E) 



CJE) , JAE) , 

^-ai(r;E) + ^^i)-(a2(r;E), 



W{E) 



W{E) 



Q^{r-E) = ^Mr-,E)-§^Mr-^E) 



W{E) 



W{E) 



where 



W{E) = ME)C3{E) - J^{E)C^{E) . 
By substituting Eq. (|4.18a|) into Eq. (|3.10|) we get to 

2m /h^ 1 



(4.18a) 
(4.18b) 

(4.19) 



G(r,s;E) 



^2m/h?E 2iJi{E) 
By substituting Eq. (|4.18b|) into Eq. ( p. 131) we get to 



a2(r;E) 



W{E) ' ^ ' W{E) 

Re(E) > 0, \m{E) > , r > s 



ai{s-E) 



(4.20) 



12 



G{r,s-E) 



2m /h^ 



C,{E) J,{E) ■ 

■0-1 (r; E) - — — o-2(r; E) 



a,{s-E) 



W{Ey'' ' W{E) 

Re(E) > 0,Im(E) < 0, r > s. (4.21) 



Since 



ai{s\E) = ai{s;E) 



(4.22) 



Eq. ( [4. 201 ) can be written as 
G(r, s; E) 



C,{E) 



, ME) 



(Ti(r; E)ai{s; E) + ^-^a2{r; E)ai{s; E) 



^2m/h?E 2iJ^{E) [ W{E) ' ' ' ' W{E) 

Re(E) > 0, Im(E) > , r > s , (4.23) 

whereas Eq. ([4.21|) can be written as 



G{r,s-E) 



2m If? 



^2m/h^E 2iJs{E) 



cm 

W{E) 



a^{r; E)ai{s; E) - 



^ ME) 



a2{r; E)a,{s; E) 



W{E) 

Re{E) > 0,lm{E) < , r > s . (4.24) 



By comparing ( |4.2|) to ( |4.23| ) we get to 

( 2m/lv' 1 -C4(E) PI 
W2ni/h^ E 2i ME)W{E) 
2mlh^ 1 1 

By comparing ( |4.2| ) to ( |4.24| ) we get to 



Re(E) > , Im(E) > 



(4.25) 



2m/ 1 C3(£^) 







nAE) 



^2m/h^E'^iME)W[E) 
2m/h^ 1 1 g 



Re{E) > , lm{E) < . 



(4.26) 



From Eqs. ( [4.25| ) and ( [4.26[ ) we can see that the measures pi2, P21 and P22 in Theorem 3 are 
zero and that the measure pu is given by 



pii{{Ei,E2)) = lim hm 



1 



E2-5 



^2 1 2m/h^ 



which leads to 



p{E) = pu{E) 



E, 47r ^2m/;i2 ^ J^{E)jm 
1 2m /tf 1 



[e+^{E - ie) - 0+ (^ + ie)] dE 
^ dE. 



47r ^2m/t9E \ME)\' 



, Eg (0,00) 



(4.27) 



(4.28) 



The function Of^{E) has a branch cut along (0, 00), and therefore (0, 00) is included in 
Sp(if). Since Sp(if) is a closed set, we have that 



Sp(iJ) = {Ei,...,E^}U[0,oo) 



(4.29) 



Therefore, the spectrum of H has a discrete part At = {-Ei, . . . , Ej^} and a continuous part 
Ac = [0,00). 
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5. DIRECT INTEGRAL DECOMPOSITION 



The fourth step toward the construction of the RHS is to compute the unitary operator 
that diagonahzes the Hamiltonian and the direct integral decomposition induced by that 
unitary operator. In order to compute them, we apply Theorem 2 of Section ^ to the 
discrete and continuous spectrum separately, since they can be treated independently. In 
the end, the independence of the discrete and continuous spectra will lead to the splitting 
of the Hilbert space into a direct sum of a Hilbert space associated to the discrete spectrum 
and another Hilbert space associated to the continuous spectrum. 

We first apply Theorem 2 to the discrete part of the spectrum. By that theorem, there 
is a unitary operator f/b from Hh '■= E(Ab)//^([0, oo), dr) onto Tib '■= defined by 

Uh ■ Hh I — > Hh 

Mr) ^ f/b/b = {(01, /b), . . . , /b)} . (5.1) 

In this equation, is the Hilbert space spanned by the bound states 0i, . . . , 0Ar, Tib is the 
space of iV-tuples of the form {{(pn, fh)}n=i (hence Tib is isomorphic to C^), /b = E(Ab)/, 
i.e., /b is the component of / along the space Tib, 

N 

Mr) = Y,{<P^,f)Ur), (5.2) 

n=l 

and (0„, /b) is the scalar product of the n-th bound state (pn{r) of Eq. ([4.11| ) with /b. 



njb)= drMr)Mr), n = l,...,iV. (5.3) 



On the space Tib, the Hamiltonian acts as the following N x N diagonal matrix: 



/El .. \ 
^2 .. 

V .. Er, J 



(5.4) 



where H\y is the restriction of H to Tib. 

Next, we apply Theorem 2 to the continuous part of the spectrum. By that theorem, 
there is a unitary map Uc from Tic := E(Ac)L^([0, oo), dr) onto L"^ {{0, oo), p{E)dE) defined 
by 

f/c : Tic I — > L\{0, oo), p{E)dE) 



/c(r) ^ ME) = {UJ,){E) = / rfr/c(r)x(r; E) , (5.5) 

where /c = E(Ac)/, i.e., /c is the component of / along the Hilbert space Tic. The space 
Tic is the Hilbert space that corresponds to the continuous part of the spectrum — hence the 
subscript c. 
If we write 
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viH) = n^,®v{nc) (5.6) 

and denote the restriction of H to T>{Hc) by H^, then the operator Uc of Eq. (^) provides a 



p-diagonahzation of He- If we seek a 5-normahzation we just have to define the following 
eigenfunction: 

<l){r;E):=^/^)x{r;E), (5.7) 

which is the eigensolution of the differential operator h that is 5-normalized, and the following 
unitary operator: 

Uc : He I — > He 



fe^ fc{E) = {UJe){E)= / rfr/,(r)0(r;E), (5.8) 

^0 

where Tic = L'^{[0,oo),dE). We note that the normalization of the bound states ( fl.ll 



IS 



different to the (5-normalization of (|5.7|) 



The inverses of the operators f/b and Uc are provided by the following theorem (cf. The- 
orem XIII. 5. 14 of Ref. 



Theorem 4 (Weyl-Kodaira) Let H, A, {pij}, etc., be as in Theorem 2. Let Eq 
and El be the end points of A. Then the inverse of the isometric isomorphism U of 
E(A)L^([0, oo), dr) onto -L^(A, {pij}) is given by the formula 



Ml 



(f/-iF)(r)= hm hm / | ^ F,(E)a,(r; E)p,,(rfi?) ) , (5.9) 



where F = [Fi, F2] G L'^{A, {pij}), the limit existing in the topology of ^^^([0, 00), dr). 
According to Theorem 4, the inverse of U^, is given by 

N 



iU^'l){r) = Y,i^n, h)Mr) , (5.10) 



n=l 



where /b = {(0i, /b), • • • , (0Af, /b)} is an A^-tuple of complex numbers. The operator U^ ^ 
transforms from Hh onto Hh- 

According to Theorem 4, the inverse of Uc is given by 

POO 

Ur) = {U-^l){r) = dEfc{E)(^{r-E), fdE) E {[0, 00), dE) . (5.11) 

Jo 

The operator transforms from He onto Tic- 
Using the operators U^ and Uc we can construct an operator U defined as 



U ■.H^H = ®^^{H{En)® / H{E)dE 

Jo 

f^Uf:= [f/b/, Ucf] = f/b, {/c(i5;)}l , (5.12) 
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where 7i is realized by Tib ® Tic, and Ti is realized by Tib ® Tic- The Hilbert space Ti(i?„), 
which is associated to each energy En in the discrete spectrum, and the Hilbert space HiE), 
which is associated to each energy E in the continuous spectrum, are realized by the Hilbert 
space of complex numbers C. On Ti, the operator H acts as the multiplication operator. 



N 



Hf ^UHf = Y, En{<Pn, /b) + {E%{E)} , / e V{H) . (5.13) 



n=l 

The scalar product on Ti can be written as 

N 



= V(/b,0n)(0n,^7b)+ r dE%{E)UE), (5.14) 

where /b = {(0i, /b), • • • , (07v, /b)} and ^b = {(0i, fl-b), • • • , (07v, S-b)}- 

Therefore, U induces a direct integral decomposition of the Hilbert space associated to the 
Hamiltonian H (see p7|,^). In particular, we can write the eigenfunction expansions of any 
element of the Hilbert spaces -Z^^([0, cxd), dr) and lv^([0, cxd), dE) in terms of the eigensolutions 
(pnif) and 0(r; i?) of /i, 

/(r) = /b(r) + /c(r) = V(0., /)0„(r) + / /(E)0(r; E) , (5.15) 

n=l 



/^/b + /c(i?) = {(0„,/)}li+ / t/r/(r)0(r;E). (5.16) 







6. RIGGED HILBERT SPACE 



As explained in the Introduction, the Hilbert space framework (in particular, the direct 
integral decomposition) is not sufficient for the purposes of Quantum Mechanics: an exten- 
sion of the Hilbert space to the RHS is needed. To construct the RHS, we shall consider the 
discrete and the continuous spectra separately. 

6.1. Construction of the Rigged Hilbert Space 

First, we construct the RHS associated to the discrete spectrum. In this case, the Hilbert 
space Tib is finite dimensional. Therefore, all the difficulties exposed in the Introduction do 
not arise. The "Gelfand triplet" that corresponds to the discrete spectrum consists of three 
copies of the Hilbert space Tib, 

*b = Tib = $b" • (6-1) 
The eigenket \En) associated to any energy E^ in the discrete spectrum is defined by 
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\En) : $b ^ C 

POO 

I — ^ {Vh\En) := (<y5b, 0n) = / dr ipb{r)(f)n{r) . (6.2) 

Clearly, the function \En) is an antilinear functional over $b, i-e., \En) G . Besides being 
an eigenvector of the Hamiltonian in the usual sense, \En) is also an eigenvector of i^b in 
the RHS sense of Eq. (|r|). 



(ip^lH^ |E„) = E„(¥^b|^n) , V(/.b e $b • (6.3) 

Although it is not necessary to use the RHS formalism when dealing with the discrete 
spectrum, we have constructed the "RHS" ( |6.1|) and the eigenket (|6l^) in order to make a 
parallel to the continuous case. 

Now, we construct the RHS associated to the continuous spectrum. The first step is to 
make all the powers of the Hamiltonian well-defined. In order to do so, we construct the 
maximal invariant subspace Vc of the operator He, 

oo 

V, := fl V{H:) . (6.4) 

n=0 

The space is the largest subspace of V{Hc) that remains stable under the action of the 
Hamiltonian and of its powers. Recalling that h denotes the formal differential operator 
( ^^ , it is easy to check that 

Dc = Wc e He I /i>c(r) G He, /i>c(0) = 0, (/.(")(«) = ^^^'Kb) = 0, n = 0, 1, 2, . . . ; 

^c(r)eC°°([0,oo))}. (6.5) 

The second step is to find the subspace C Vc on which the eigenkets \E) of Eq. ( |1.5b| ) 
act as antilinear functionals. Those eigenkets are to be defined as integral operators whose 
kernel is the eigenfunction 0(r; E), 

\E) : $c ^ C 



cp, ^ {ip,\E) := dr ip,{r)(l){r; E) = {U,ip,)iE) . (6.6) 

^0 

A necessary condition for \E) to be well-defined is that $c be included in Vc- This con- 
dition, however, is not sufficient to obtain a continuous functional. In order to make the 
eigenfunctional \E) continuous, we have to impose further restrictions on the behavior at 
infinity of the elements of V^, 

rfr |(r + l)"(/i + l)'"(^c(r)|^ < oo, m = 0, 1, 2, . . . (6.7) 



The imposition of (pTZD upon the space Vc yields 



POO 

{ifcEV^l / rfr |(r + l)"(/i+l)>e(r)|' < cx), n, m = 0, 1, 2, . . .} . (6.8) 
Jo 
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On $c, we define the family of norms 



lbc||n,m:=yy rfr |(r + l)"(/i + l)-^e(r)|', n, m = 0, 1, 2, . . . (6.9) 

These norms can be used to define a countably normed topology on $c (see [Q), 

^c,a-^^^^c iff llv^c.a - V^clU.m ^0, n, m = 0, 1, 2, . . . (6.10) 



a— >oo a— >oo 



Once we have constructed the space we can construct its topological dual as 
the space of r^^-continuous antilinear functionals on $c (see 0]) and therewith the RHS 
corresponding to the continuous part of the spectrum, 

*c C K C . (6.11) 



In order to show that the RHS (|6.11| ) is what we are looking for, we need to prove the 
following proposition (the proof can be found in Appendix P): 

Proposition 1 The triplet of spaces (|6.11|) is a Rigged Hilbert Space, and it satisfies 
all the requirements demanded in the Introduction. More specifically. 



(i) The quantities ( p.9|) fulfill the conditions to be a norm. 

(ii) The space $c is stable under the action of iJc, and is r^^-continuous. 



(iii) The ket \E) of Eq. (|6.6|) is a well-defined antilinear functional on i-6., \E) belongs 

to 

(iv) The ket \E) is a generalized eigenvector of ifc, 

H:\E) = E\E), (6.12) 

i.e., 

{^c\H:\E) = {Hl^c\E) = E{ip,\E) , V^e e (6.13) 



(This proposition is a generalization to infinite dimensional spaces of the results that 
hold in finite dimensional spaces.) 

We now combine the RHSs ( |6.1D and ( |6.11D corresponding to the discrete and continuous 
spectrum into the RHS of the square well-barrier potential, 

$c7i:c$\ (6.14) 

where 

$:=$b©^c; 'H:='H^®n,] $^:=$^©$^ (6.15) 

On the space all the expectation values of the Hamiltonian and all the algebraic operations 
involving H are well-defined, and the eigenvalue equations ( |1.5|) hold. The kets |i?„) and \E) 
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satisfy the normalization ( |1.6|) . Since the spaces $b and $c are orthogonal to each other, 
we also have that 



{^\En) = (<^b|K) ; {v\E) = {ip,\E) . (6.16) 

As we shall see in the next section, the kets \En) and \E) form a complete basis system. 

Before finishing this section, we would like to remark that 0(r; E) is not the same object 
as the ket \E) of Eq. ( |6.6| ). The function 0(r; E) is an eigenfunction of the formal differential 
operator h, whereas the ket \E) is a generalized eigenvector of the Hamiltonian. The eigenket 
\E) is defined by the eigenfunction 0(r; E) and by the space of test functions #c on which 
it acts as an antilinear functional. 



6.2. Dirac Basis Vector Expansion 

In the Introduction, we mentioned that in Quantum Mechanics it is assumed that the 
eigenkets of the Hamiltonian form a complete basis system that can be used to expand any 
physical wave function ip as in Eq. ( |1.7| ). This expansion is derived in the present section. 
That derivation consists of the restriction of the eigenfunction expansion ( ^.15| ) to the space 

If we denote {r\(f) = ip{r), {r\En) = 0n('"), and {r\E) = (f){r;E), and if we define the 
action of the bra {E\ on the wave function ip by {E\(p) := {ip\E), then Eq. ( [5.15[ ) can be 
written as 

{r\ip) = J2{r\En){E^\^)+ dE {r\E){E\ip) , ^G$. (6.17) 

n=0 

This is the Dirac basis vector expansion of the square well-barrier potential. Although the 
eigenfunction expansion ( |5.15|) is valid for every element of the Hilbert space, Dirac basis 
vector expansion ( |6.17| ) is only valid for functions because only those functions fulfill 

both 



cp^iE) = {ip\E) (6.18) 

and 

{(p\H''\E) = {W^\E) = E{y^\E) . (6.19) 

For the sake of completeness, we include the Nuclear Spectral Theorem [^, which is 
usually referred to as the mathematical justification of the heuristic Dirac basis vector ex- 
pansion [|l]-^. Instead of using the general proof of Ref. 0], we prove that theorem using 
the machinery of the Sturm-Liouville theory (our proof can be found in Appendix 0) . 

Proposition 2 (Nuclear Spectral Theorem) Let 

$ C L2([0,oo),dr) C (6.20) 

be the RHS of the square well-barrier Hamiltonian H such that $ remains invariant under 
H and if is a r^-continuous operator on Then, for each energy in the spectrum of H 
there is a generalized eigenvector such that 



19 



H\En) = En\En) , En e {El, . . . , ^iv} , (6.21a) 
H''\E) = E\E), Eg[0,cx)), (6.21b) 



and such that 

N 



POO 

((^,^) = V(¥.|EO(i?n|^)+ / dE{^\E){E\^), Vy.,^e$, (6.22) 



n=l 

and 

N 



" /-oo 
(^,if™V^)=5^Er(¥'|i^n)(i?n|V^)+ / dEE^{cp\E){E\^l,), V^,^e$^ 



m = 1,2, 



(6.23) 

Thus this theorem allows us to write the scalar product ( |6.22| ) of any two functions ip 



of $ and the matrix elements ( |6.23|) in terms of the action of the kets \En) and \E) on (f^tl^. 



6.3. Energy and Momentum Representations of the Rigged Hilbert Space 

In this section, we construct the energy representation of the RHS (|6.14| ) by applying 
to it the unitary operator U of Eq. ( |5.12|) . As done throughout this paper, we study the 
discrete and the continuous case separately. 

In the discrete case ( |6.1|) , we have already seen that the energy representation of Hh is 
C^. Hence the energy representation of the RHS ( p-l\) consists of three copies of the space 
C^. The image of the eigenket \En) under is the A^-tuple that has zeros everywhere but 
at the n-th position, {Sin}fLi. In the energy representation, the Hamiltonian i^b acts as the 
diagonal matrix ( |5.4| ). 

In the continuous case (|6.11 ), we have already shown that in the energy representation 



the Hamiltonian He acts as the multiplication operator. The energy representation of the 
space $c is defined as 

*c := f/c$c . (6.24) 

The space $c is a linear subspace of L'^{[0, oo), dE). In order to endow $c with a topology 
rg^, we carry the topology on $c into 

■■= UcT^. ■ (6.25) 

With this topology, the space $c is a linear topological space. If we denote the dual space 
of $c by , then we have 

U,^^: = {U,^,r = . (6.26) 
From the following chain of equalities: 
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(^,|f/,X|i?) = (f/-Vc|i?) 



dr ipc{r)(j){r; E) 



= Vc{E) , (6.27) 

it follows that the energy representation of the Dirac ket \E), that we denote by \E), is the 
antilinear Schwartz delta functional, i.e., \E) is the antilinear functional that associates to 
each function (p^ the complex conjugate of its value at the point E. 

In order to summarize all these results, it is very helpful to show the different realizations 
of the RHS through the following diagram: 

H- if $b © $c C T^b © C © \En) © 1^) 

t/b i © i f/c f/b i © i f/c f/b i © i (6.28) 

^ $b © C Tib © C $b © ^c^ \En) © 1^) 



On the top line of the diagram ( 5.28 ), we have the RHS, the Hamiltonian, the wave functions 
and the Dirac kets in the position representation. On the bottom line, we have their energy 
representation counterparts. 

Since every function we have computed (eigenfunctions of the Schrodinger differential 
operator. Green function, etc.) depends on the square root of the energy rather than on the 
energy itself, the momentum is a more convenient variable, especially for the continuous 
case. The momentum k is defined by 



1 2m 

k:=\l — E. (6.29) 



In terms of /c, the 5-normalized eigensolution of the differential operator ( p. 71) reads 

{r\k) = [2^jmMk)]-''^ x{r; k) . (6.30) 

This eigensolution can be used to obtain the unitary operator Vc that transforms from the 
position into the momentum representation. 



fc{k) = {VJ,){k) = drf,{r){r\k), f,eH,. (6.31) 



(The operator Vt, corresponding to the discrete case is f/b-) On the ^-representation, the 
(continuous part of the) Hamiltonian acts as multiplication by k"^. To each k G [0, oo), there 
corresponds an eigenket \k) that acts on $c as the following integral operator: 



{ip,\k) := / dr {ipc\r){r\k) = (K<^c)(A;) , <^c G $c ■ (6.32) 

^0 

This eigenket is a generalized eigenvector of the Hamiltonian corresponding to the eigenvalue 
k'^. In this way, we can construct the ^-representation of all the objects of the diagram ( |6.28D . 



If we express our functions in terms of k, there is no need to distinguish different regions in 
the ii^-plane, because in the fc-plane the resolvent and the S*- matrix do not have a cut. There 
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is no need to work with a particular branch of the square root either, because Eq. (|6.29|) 
provides a Riemann surface in a natural way. For instance, we can define the Jost function 



J{k) := J+{k) , keC, (6.33) 

the eigensolution 

/(r;A;) :=e+(r;fc), keC, (6.34) 

and the Green function 

' i\ 2m/ff x{r<;k)f{ry,k) 
G(r, r;fc):= , — , keL, (6.35) 

where r<, r> refer to the smaller and to the bigger of r and r', respectively. Then, as 
k approaches the real positive k axis, or as E approaches the right-hand cut from above, 
i7(fc), f{r;k) and G{r,r';k) become J'+{k), 9+(r; fc) and G~^{r,r'; k), respectively; as k 
approaches the negative real axis (from above), or E the right-hand cut from below, J^{k), 
f{r;k) and G{r,r';k) become J_{k), 6_(r;A;) and G~ {r,r'; k), respectively |2T|. Also, the 



values that the functions of Eqs. ( |6.33|) - (|6.35|) take on the second sheet of the Riemann 
surface are already specified (the Sturm-Liouville theory deals only with energies in the first 
sheet). This is particularly useful when studying resonance energies, which lie on the second 
sheet [El-El. 



7. CONCLUSIONS 

We have shown that the solutions of the Schrodinger equation corresponding to the 
square well-barrier potential fall in the Rigged Hilbert Space 

$b © *c C T^b © C © . (7.1) 

The spectrum of the Hamiltonian has a discrete part {Ei, . . . , E^} and a continuous part 
[0,cxd). For each energy in the spectrum of H, we have constructed an eigenket that is 
an eigenvector of H. If the energy En belongs to the discrete spectrum, its corresponding 
eigenket |-E„), which is given by Eq. ( |6.2| ), is an element of If the energy E belongs to 
the continuous spectrum, its corresponding eigenket \E), which is given by Eq ( |6.6| ), is an 
element of . Thus the RHS treats the discrete and the continuous spectra independently 
and on the same footing. Each element (/? of © can be expanded by the eigenkets 
of the Hamiltonian as in Eq. ( |6.17 ). The elements of $b © are the only square normal- 



izable functions for which such an expansion is possible. The expectation values and the 
uncertainties of the Hamiltonian are well defined quantities in each element (f of $b © ^c- 

Therefore, the RHS ( |7. 1|) contains all the physically meaningful solutions of the 
Schrodinger equation — the physical wave functions are included in $b©^c, whereas the mo- 
noenergetic solutions of the time independent Schrodinger equation are included in ©$^. 

Although illustrated within the example of the square well-barrier potential, these con- 
clusions hold in general: the solutions of the Schrodinger equation fall in a RHS rather than 
just in a Hilbert space. 
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APPENDIX A: LIST OF AUXILIARY FUNCTIONS 



If we define 



, 2m ^ ~ 
k:=\l-—E- Q,:- 



2m 



2m 



{E - V2) . 



then the coefficients in Eq. (p.7|) are given by 



J^{E) = J,{k) 



ME) ^ J,{k) = ^e^^"^ 



JsiE) = Js{k) 
ME) = Mk) 



4 
1 

r 



-kb 



1 + Si 1 e«i" 
Q2 

1 - Si 1 e«^" 
Q2 



1 _ Si 1 e-'^^" 
Q2 

1 + Si ) e-^^^ 
Q2J 

Q2 



1 + ^] e^-''Mk) + 1 - :^ e-'^'^'Mk) 



k 



k 



\ kb 

2 



1 - S j e^^'Mk) + I 1 + y 1 e-^^'Mk) 



Q2 



and the coefficients in Eq. ( ^.9| ) by 
1 r k 



AsiE) = Asik) 
ME) ^ Mk) 
ME)=Mk) 



-Q2b I X _ ^ , g 



■kb 



Q2 



l-e^'^' I 1 + i 1 e 

2 V Q^^ 



kb 



1 



1 + S^ ) e^-^Asik) + ( 1 - £ 1 e^'^-'^A.ik) 



Qi 



Qi 



A2{E) ~ A(A:) = ^e'^^" 



If we define 



then the functions JiE) of Eq. (|3.11|) are given by 



(Al) 



(A2a) 
(A2b) 
(A2c) 
(A2d) 



(A3a) 
(A3b) 
(A3c) 
(A3d) 



(A4) 
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2 V 



J2{E) = J2{k) = -e^^^'^ f sin(gia) - ^ cos(Qia; 



-ikb 



1 + 



ME) = Mk) = ^e'^"^ 



9i 
k 



02 
k 



and the functions A+{E) of Eq. ( |^ ) by 



At{E)^At{k) = \e-^^^' 



A; 



At{E)=At{k) = -e-^Qi'^ 



1 + e^Q^'^^+(fc) + f 1 - ^ ) e-''^'''At{k) 



Qi 



Q2 



The coefficients in Eq. ( ^.141 ) are given by 



k 



-ikb 



A,{E)^A^{k) = -e'^^^' 



A^{E)^A^{k) = -e-^^^^ 
A^{E)^A^{k) = ^e^<^^^ 



JQ2a A-n.\ , f, , Q2 



The coefficients in Eq. ( [4.4aD are given by 



(A5a) 
(A5b) 
(A5c) 
(A5d) 

(A6a) 
(A6b) 
(A6c) 
(A6d) 

(A7a) 
(A7b) 
(A7c) 
(A7d) 



B,{E) ^ B,(k) = ^e^-' fl - ^ 1 e"' 



B,{E)=B,{k) = ^e-'^-- 



1 + 2^ 1 e'^-'^Bsik) + I 1 



Q2 
Qi 



-Q2a 



B,{k) 



Q 
Qi 



2 1 ^Q2a 



e'^'^Bsik) + 1 + 



Q2 
Qi 



-^^"Biik) 



(A8a) 
(A8b) 
(A8c) 
(A8d) 



The functions C{E) of Eq. ([4.1 7b| ) are given by 
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C,{E)=C,ik) 

C2{E)=C2{k) 

C,{E)^C,{k) 
C,{E) = C,{k) 



ie-^Q^'^ (cosiQ^a)-■^sm{Q^a: 
2 V 

1 



-ikb 



cos{Qia) + sm{Qia 



iQ2 

1 + e'Q^'Cik) +fl-^) e-'^^'^ik) 



ikb 



k 



l + ^)e-^Q%{k) 



(A9a) 
(A9b) 
(A9c) 
(A9d) 



n,m — llV-'clln.m ~l~ ||V^c||n,m ) 
II '-'^V-'c II n,m I'-'^l II V-'c II n,m ) 
II V-'c II n,m ^ , 

If Ilv^c||n,m = 0, then v^c = . 



APPENDIX B: PROOFS OF PROPOSITIONS 1 AND 2 

Proof of Proposition 1 . 

(i) It is very easy to show that the quantities fulfill the conditions to be a norm, 

(Bla) 
(Bib) 
(Blc) 
(Bid) 

The only condition that is somewhat difficult to prove is (Pld|) : if | |(y3c 1 1 n,m = 0, then 

{l + rr{h+ir^,{r)=Q, (B2) 

which yields 

(/i+l)>e(r) = 0. (B3) 

If m = 0, then Eq. ( p3|) implies ^d''^) =0. If m = 1, then Eq. ( P3|) implies that —1 is an 
eigenvalue of He whose corresponding eigenvector is (fc- Since there is no discrete eigenvalue 
in the continuous part of the spectrum, ip^ must be the zero vector. The proof for m > 1 is 
similar. 

(ii) In order to see that is r$^-continuous, we just have to realize that 

1 1 -f^cV^c 1 1 n,m = \\{.Hc + I)'Pc — Vc\\n,m 

— II (-^c ~l~ -^)V-'c II n,m ~l~ II V-'c II n,m, 

= ll'/'clln.m+l + llV^clU.m • (B4) 

The stability of $c under the action of He also follows from Eq. 



(iii) From the definition (|6.6|) , it is pretty easy to see that \E) is an antilinear functional. 
In order to show that \E) is continuous, we define 



M{E) := sup |0(r;E)| . 

re[0,oo) 



(B5) 
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Since 



\{v>c\E)\ 



dr (y9c(r)0(r; E) 



<M{E) / dr\^^{r) 



M{E) r dr-^{l + r)\^,{r)\ 
Jo 1 + ^ 



<M{E) / dr 



A^(^)ll<^c||i,o, 



1/2 



(ir 1(1 + r)ipc{r)[ 



1/2 



(B6) 



the functional \E) is continuous when $c is endowed with the topology. 



(iv) In order to prove that \E) is a generalized eigenvector of He, we make use of the 
conditions ( |6.5|) and ( |6.7| ) satisfied by the elements of 

{cf,\H^\E) = {Hl^,\E) 

d^ 



dr 



2m dr'^ 



+ Vir) if,{r)(f){r;E) 



2m 



dr 



[r-E) 



+ 



+ / dnpcir) 
Jo 

E{^,\E). 





2m dr"^ 



2m 



rf0(r; E) 
dr 



J 



+ V{r) (P{r-E) 



(B7) 



Proof of Proposition 2. 

Let if and ■?/' be in Since f/ of Eq. (|5.12| ) is unitary, 
The scalar product of and ^/'b can be written as 

N 



(^b,V'b)Wb = '^{(p\En){E„ 



n=l 



(B8) 



(B9) 



The wave functions ipc and i/Jc are in particular elements of L'^{[0,oo),dE). Hence their 
scalar product is well-defined, 



(^c,V'c)we= / dE^,{E)ME). (BIO) 

Since and ipc belong to the action of each eigenket \E) on them is well-defined, 

(Blla) 
(Bllb) 



{^,\E) = (p,{E) , 



By plugging Eq. ( pTID into Eq. (|BT0D and Eqs. ( |BT0|) and (^ into Eq. 0, we get to 
Eq. ( 16.221) . The proof of (|6.23|) follows the same pattern. 
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